Spatial inhomogeneity and strong correlation physics: a dynamical mean field study of 
a model Mott-insulator/band-insulator heterostructure 
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We use the dynamical mean field method to investigate electronic properties of heterostructures 
in which finite number of Mott-insulator layers are embedded in a spatially infinite band-insulator. 
The evolution of the correlation effects with the number of Mott insulating layers and with position 
in the heterostructure is determined, and the optical conductivity is computed. It is shown that 
the heterostructures are generally metallic, with moderately renormalized bands of quasiparticles 
appearing at the interface between the correlated and uncorrelated regions. 
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An exciting new direction in materials science is 
the fabrication and study of heterostructure involv- 
ing "correlated electron " materials such as Mott in- 
sulators, high-temperature superconductors, and novel 
magnets 0, The issues raised by these heterostruc- 
tures, especially the evolution with position of proper- 
ties from correlated to uncorrelated, is of fundamen- 
tal physical interest and would be crucial for prospec- 
tive devices based on correlated electron compounds. 
Many interesting systems have been fabricated, includ- 
ing modulation-doped high-T^ superconductors 0, 0, 
Mott-insulator/band-insulator heterostructures 0, and 
a variety of combinations of magnetic transition-metal 
oxides 0, 0, but there has been relatively little the- 
oretical study of the heterostructure-induced changes in 
many-body physics. The theoretical problem is difficult 
because it requires methods which can deal both with 
spatial inhomogeneity and strong correlation physics. 

Despite the difficulties, several interesting works have 
appeared. Fang, Soloviev and TerakuraQ used bulk 
band structure calculation to gain insight into the effects 
of strain fields induced by lattice mismatch in a het- 
erostructure. Matzdorf and co-workers used GGA band 
theory methods to study the surface electronic and lat- 
tice structure of Sr2Ru04, predicting a lattice distortion 
which is observed and surface ferromagnetism which is 
apparently not observed Potthoff and Nolting used 
dynamical-mean-field methods to study the consequences 
of the lower coordination at a surface 0, 0| . 

All of these papers, however, treated situations in 
which the electronic density remained at the bulk value, 
and the new physics arose from structural differences. 
A crucial feature of heterostructures is an inhomoge- 
neous electron density caused by a spreading of charge 
across the interfaces which define the system. Recently 
[13I H^ . we used realistic multi-orbital interaction pa- 
rameters and a density-functional-theory-derived tight- 
binding band structure to model ground state proper- 
ties of the LaTiOa/SrTiOs heterostructure fabricated by 
Ohtomo et al.^. While this study captured important 
aspects of the density inhomogeneity, it did not address 



the dynamical properties of correlated heterostructures. 
Further, this study employed the Hartree-Fock approx- 
imation which is known to be an inadequate represen- 
tation of strongly correlated materials, and in particular 
does not include the physics associated with proximity 
to the Mott insulating state. 

In this paper we use the dynamical mean field method 
which provides a much better representation of the 
electronic dynamics associated with strong correlations, 
to study the correlated electron properties of a simple 
Hubbard-model heterostructure inspired by-but not a 
fully realistic representation of-the systems studied in 
We present results for observables including photoe- 
mission spectra, optical conductivity, charge density, and 
highlight similarities and differences to previous work. 

We study a single orbital model with basic Hamilto- 
nian Hhub ~ Hband + Hint + Hcoui with 
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Here the sites i form a simple cubic lattice of lattice con- 
stant a, so electronic positions Ri — a{ni,mi,li). We 
include both an on-site (J7) and long ranged Coulomb 
interaction: the screening field from the latter is im- 
portant for the electron density profile. We associate 
the electronic sites with the B-sites of an ABO3 per- 
ovskite lattice, and define the heterostructure by coun- 
terions of charge -1-1 placed on a subset of the A-sites. 
Here we study an n-layer [001] heterostructure defined 
by n planes of -|-1 counterions placed at positions = 
a{nj + 1/2, nij + 1/2, Ij + 1/2), with —00 < nj, nij < 00 
and the Ij runing over n adjacent values. The resulting 
potential is 
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Charge neutrality requires that the areal density of elec- 
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trons is n. A dimensionless measure of the strength of 
the Coulomb interaction is Ec = /{eat); we choose pa- 
rameters somewhat arbitrarily so that = 0.8 (this cor- 
responds to t ~ 0.3 eV and length a ~ 4 A and e — 15, 
which describe the system studied in Q). We found that 
the charge profile did not depend in an important way 
on e for 5 < e < 25. 

The basic object of our study is the electron Green 
function, which for the [001] heterostructure may be writ- 
ten 

G{z,z' ,k\\\Uj) = [uj+ii-Hband-Hcoui-'^{z,z' ,k\\\Uj)]~'^. 

(4) 

We approximate the self-energy operator as the sum of 
a Hartree term arising from the long-ranged part of the 
Coulomb interaction 



£\Ri — Rn 
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and a dynamical part E^) arising from local fluctuations . 
Following the usual assumptions of dynamical-mean-field 
theory |l5| as generalized to inhomogeneous situations by 
Schwieger et. al. we assume 



(6) 



The layer (z)-dependent dynamical self energy is de- 
termined from the solution of a quantum impurity model 
with mean-field function fixed by the self-consistency 
condition 



G""P(z,cj) 
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One must solve a separate impurity model for each layer, 
but the self consistency condition [cf. Eq. Q] implies 
the solutions are coupled. It is also necessary to self- 
consistently calculate the charge density via n{z) = 
-2 / ^/^ImG™P(2,w) with / the Fermi distribution 
function. The numerics are time consuming, and it is 
therefore necessary to adopt a computationally inexpen- 
sive method for solving the quantum impurity models. 
We use the two-site method of Potthoff jla|, which re- 
produces remarkably accurately the scaling of the quasi- 
particle weight and lower Hubbard band near the Mott 
transition. We have also verified 0| that the two site 
method reproduces within a few percent the T — mag- 
netic phase diagram found by Ulmke in ^ model with 
an unusual low-energy density of states peak. 

Fig. n shows the layer-resolved spectral function 

A{z,z;uj) — — i / jj:^lmG{z, z,k\i;Lj + iO+) for a 10- 
layer heterostructure with U — 16t (about 25% greater 
than the critical value which drives a Mott transition 
in a bulk system described by Hhut with n = 1. The 
spectral functions are in principle measureable in pho- 
toemission or scanning tunneling microscopy. Outside 




FIG. 1: Layer-resolved spectral function calculated for 10- 
layer heterostructure for U — 16t, e = 15. The heterostruc- 
ture is defined by -1-1 charges placed at z — ±0.5, ±1.5, ...±4.5 
so the electronic (B) sites are at integer values of z. 



the heterostructure (z > 6), the spectral function is es- 
sentially identical in form to that of the free tight-binding 
model Hhand- The electron density is negligible, as can be 
seen from the fact that almost all of the spectral function 
lies above the chemical potential. As one approaches the 
heterostructure (z = 6), the spectral function begins to 
broaden. Inside it (z < 5) weight around w = begins to 
decrease and the characteristic strong correlations struc- 
ture of lower and upper Hubbard bands with a central 
quasiparticle peak begins to form. The sharp separation 
between these features is an artifact of the 2-site DMFT 
[as is, we suspect, the shift in energy of the upper (empty 
state) Hubbard band for z = 4,5]. Experience with bulk 
calculations suggests that the existence of three features 
and the weight in the quasiparticle region are reliable. 
Towards the center of the heterostructure, the weight 
in the quasiparticle band becomes very small, indicating 
nearly insulating behavior. For very thick heterostruc- 
tures, we find the weight approaches exponentially. 

The behavior shown in Fig. ^ is driven by the varia- 
tion in density caused by leakage of electrons out of the 
heterostructure region. Fig. |21 shows the numerical re- 
sults for the charge-density distribution n(z) for the het- 
erostructure whose photoemission spectra are shown in 
Fig. n One sees that in the center of the heterostructure 
(z = 0) the charge density is approximately 1 per site, 
and that there exists an edge region, of about three-unit- 
cell width, over which the density drops from ~ 1 to ~ 0. 
The over-all charge profile is determined mainly by the 



3 



§ 0.6 
=3 



15 0.4 




FIG. 2: Total charge density (open squares) and charge 
density from coherent part near Fermi level (open circles). 
For comparison, total charge density calculated by apply- 
ing Hartree Fock approximation to the Hamiltonian shown 
as filled symbols. Parameters as in Fig. Q 



self consistent screening of the Coulomb fields which de- 
fine the heterostructure, and is only very weakly affected 
by the details of the strong on-site correlations (although 
the fact that the correlations constrain n < 1 is obviously 
important). To show this, we have used the Hartree-Fock 
approximation to recalculate the charge profile: the re- 
sults are shown as filled circles in Fig. [3 and are seen to 
be almost identical to the DMFT results. 

The existence of an approximately three unit cell wide 
edge region where the density deviates significantly from 
the values n = and n — 1 characteristic of the two sys- 
tems in bulk form implies that only relatively thick het- 
erostructures (n > 6) will display 'insulating' behavior 
in their central layers, and suggests that the edge regions 
sustain quasiparticle subbands which give rise to metallic 
behavior. The open circles in Fig.[21show the charge den- 
sity in the 'quasiparticle bands' [obtained by integrating 
A{z,z;u;) from uj — down to the first point at which 
A{z,z;uj) = 0]. One sees that these near Fermi-surface 
states contain a small but non-negligible fraction of the 
total density, suggesting that edges should display rel- 
atively robust metallic behavior. The results represent, 
a significant correction to the Hartree-Fock calculation 
[l^ . which leads, in the edge region, to a metallic quasi- 
particle density essentially equal to the total density. 

The spectral function is determined by the layer- 
dependent dynamical self-energy TtD{z,uj). In bulk ma- 
terials one distinguishes Fermi liquid and Mott insula- 
tors by the low- frequency behavior of T,o', in a Fermi 
liquid TtD -^uj^o (1 — Z~^)uj (leading to a quasiparti- 
cle with renormalized mass) while in a Mott insulator 
T,D ^ui^o A^/w (leading to a gap in the spectrum). 
In the heterostructures we study, we find that outside 
the high density region, correlations are weak {Z « 1), 
and that as one moves to the interior of thicker het- 
erostructures, correlations increase {Z decreases). Mott 
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FIG. 3: (a) Dispersion relations of fiUed-subband quasi par- 
ticles calculated for 3 layer heterostructure with U = 16i and 
£ — 15. Solid and broken lines are for odd- and even-a, respec- 
tively, (b) Subband quasiparticle weights, (c, d) Quasiparti- 
cle wave functions for q = 1, 3, 5 and 2,4 at w = 0. Here the 
heterostructure is defined by -1-1 charges placed at z = 0, ±1 
so the electronic (B) sites are at half integer values of z. 



insulating solutions {Z = 0) are never found; instead 
J^d{z,i^) ~ [1 - Z~^{z)]uj with < Z < 1 for all layers 
z, although in the interior of thick, large U heterostruc- 
tures Z is only nonvanishing because of leakage (quan- 
tum tuneling) of quasiparticles from the edges, and goes 
exponentially towards zero. 

The nonvanishing Z indicates a Fermi liquid state with 
well defined coherent quasiparticles (thus negligible low 
frequency scattering). In the heterostructure context the 
quasiparticles form subbands, with quasiparticle energies 
Ea(k\\) and wave functions (fa{z', Ea{k\\)) which are the 
low energy eigenfunctions and eigenvalues of 



Z ^{z)EaS^,y + ^i- Hband{k\\) " Hcoul 'fa{z') 



0. 

(8) 

(We note that the 2-site DMFT method used here is be- 
lieved to give reasonable results for Z but of course ne- 
glects scattering effects. Near Ea = scattering is unim- 
portant but of course will increase at higher energies.) 

Numerical results for the coherent quasiparticles in a 
heterostructure with n = 3 and U — IQt are shown in 
Fig. 121 For these parameters we find 8 quasiparticle 
bands with non-vanishing electron density. The calcu- 
lated dispersion relations are shown in panel (a) and are 
labelled a = 1...8 in order of decreasing electron density. 
We observe that the band splittings depend on momen- 
tum because of the layer dependence of Z . The corre- 
sponding quasiparticle weights Za and real-space wave 
functions ipa{z) at w = are shown in Fig. Elb) and 
(c,d), respectively. (These quantities vary somewhat over 
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FIG. 4: Heavy lines: low frequency (quasiparticle-region) op- 
tical conductivity for n = 3, U — 16t heterostructure. Peaks 
arise from transition among subbands shown in Fig. |3Ja). 
Light lines: conductivity calculated from Hartree-Fock ap- 
proximation to same Hamiltonian. 

the band also). Za is the smallest for the a — 1 sub- 
band because its real-space wave function contains the 
largest weight ai z = ±0.5, where the charge density is 
the largest and, therefore, the correlation effect is the 
strongest [see Fig. EIc)]. Za generally increases with 
increasing a, because as a increases the wave function 
amplitudes |(y5Q,(z)| decrease in the high density regions 
(near z = 0). The anomalies observed in Za at a = 4, 6 
correspond to the increase of \(fa{z = ±0.5) | due to the 
symmetry of the wave function [see Fig.EJd)]. 

The coherent subbands may be studied by optical con- 
ductivity with electric field directed along [001]. As an 
example, the heavy line in Fig. 0| shows the quasipar- 
ticle contribution to the conductivity spectrum, calcu- 
lated for a heterostructure with n = 3 and U = 16t 
using the standard Kubo formula with optical matrix el- 
ement obtained by applying the Peierls phase ansatz to 
Hband (and t = 0.3 eV). Three main features are ev- 
ident dX Ld — 0.75t, 0.55< and 0.45i; each of these has 
contributions from two interband transitions, (1 — > 2, 
2^3), (3 ^ 4, 4 ^ 5) and (5 ^ 8, 6 ^ 7) respectively. 
The optical features are not sharp because the quasipar- 
ticle band splitting depends on . The weaker features 
at lower energies arise from transitions involving high- 
lying, only slightly occupied bands. The lighter lines in 
Fig. 21 show the optical conductivity computed using the 
Hartree-Fock approximation. We see that the spectra 
are qualitatively similar, but that the Hartree-Fock ab- 
sorption features occur at a larger energy because the Z- 
induced band narrowing is absent and are delta functions 
because in the Hartree-Fock approximation the subbands 
splittings are fc|| -independent. 

To summarize, we have presented the first dynamical 
mean field study of a 'correlated electron heterostruc- 
ture', in which the behavior is controlled by the spread- 
ing of the electronic charge out of the confinement region. 



Our results show how the electronic behavior evolves 
from the weakly correlated to the strongly correlated re- 
gions, and in particular confirms the existence of an ap- 
proximately three unit cell wide crossover region in which 
a system, insulating in bulk, can sustain metallic behav- 
ior. We found that even in the presence of very strong 
bulk correlations, the metallic edge behavior displays a 
finite (roughly factor-of-two-to-three) mass renormaliza- 
tion. We showed how the magnitude of the renormaliza- 
tion is affected by the spatial structure of the quasiparti- 
cle wave function and determined how this renormaliza- 
tion affects physical properties, in particular the optical 
conductivity. Important future directions for research is 
to examine the phase diagram, using beyond Hartree- 
Fock techniques, and to generalize the results presented 
here to more complicated and realistic cases. 
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